Quantum information processing in self-assembled crystals of cold polar molecules 
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We discuss the implementation of quantum gate operations in a self-assembled dipolar crystal 
of polar molecules. Here qubits are encoded in long-lived spin states of the molecular ground 
state and stabilized against collisions by repulsive dipole-dipole interactions. To overcome the 
single site addressability problem in this high density crystalline phase, we describe a new approach 
for implementing controlled single and two-qubit operations based on resonantly enhanced spin- 
spin interactions mediated by a localized phonon mode. This local mode is created at a specified 
lattice position with the help of an additional marker molecule such that individual qubits can 
be manipulated by using otherwise global static and microwave fields only. We present a general 
strategy for generating state and time dependent dipole moments to implement a universal set 
of gate operations for molecular qubits and we analyze the resulting gate fidelities under realistic 
conditions. Our analysis demonstrates the experimental feasibility of this approach for scalable 
quantum computing or digital quantum simulation schemes with polar molecules. 

PACS numbers: 03.67.Lx, 33.80.Ps, 61.50.-f 



I. INTRODUCTION 

Rapid progress in control of atomic, molecular and op- 
tical (AMD) systems in the recent past has established an 
experimental basis to realize quantum information pro- 
cessing (QIP) schemes with great success [iHl]. All basic 
ingredients that are necessary for QIP [5], i.e. highly 
developed state preparation, single particle control and 
measurement as well as controlled many body gate oper- 
ations have been experimentally demonstrated with cold 
atoms and ions, which take a pioneering role in this con- 
text. Prominent examples include ion trap experiments, 
from early realizations of controlled gate operations |6j 
to recently achieved complex quantum simulations [7], 
as well as the addressing of single spins [8] and a variety 
of many body gate implementations [9l [10] in systems 
of neutral atoms. In recent years, parallel to this devel- 
opment considerable theoretical and experimental efforts 
have been focused on achieving a similar level of con- 
trol over polar molecules pJTfTT] . While state prepara- 
tion and manipulation of single molecules are still part 
of ongoing efforts, techniques for the generation of dipo- 
lar gases in the quantum degenerate regime have already 
been demonstrated in the lab [TSHST]. The unique fea- 
tures of polar molecules [4 , like controllable permanent 
dipole moments and long lived rotational states with ad- 
dressability by microwave fields, offer new and interesting 
possibilities for QIP. 

The first quantum computation scheme with polar 
molecules has been proposed by Dave DeMille [22 with 
molecules trapped at separate sites of a ID optical lattice 
and qubits encoded in the rotational states. In this setup 
individual site addressability is achieved by spectrally re- 
solving distinct qubits using an electric gradient field and 
a strong and state dependent dipole-dipole interaction 
between neighboring molecules is employed to implement 



two qubit gates. To extend this basic concept beyond an 
"always-on" interaction using switchable dipole moments 
has been proposed p3j with resemblance to QIP con- 
cepts with Rydberg atoms in optical lattices [2T. Other 
ideas include hybrid quantum computing schemes with 
atom-molecule platforms p5l [26j^ or molecules coupled to 
stripline cavities [27'J30|. 

In contrast to preceding proposals here we want to fo- 
cus on QIP in self- assembled molecular dipolar crystals 
(MDC). Such (quasi-) crystals form under strong repul- 
sive dipole-dipole interactions in an ensemble of polar 
molecules confined to a plane with perpendicular dipole 
moments aligned by a strong electric field [3TH33] . Con- 
trary to Wigner crystals that rely on Coulomb interac- 
tions the formation of an MDC is favored by high densi- 
ties due to the fast decay of the dipole-dipole interaction 
and lattice spacings range down to a ~ 100 nm [3T] . 
From a perspective of QIP this would be advantageous 
as it leads to potentially stronger molecular interactions 
in comparison to previously described setups with opti- 
cal lattices. In addition, the self-assembled crystal pro- 
vides a stable defect-free periodic structure even at finite 
temperatures up to ~ /iK. Finally, the trapping poten- 
tial for a MDC only requires a strong transverse confine- 
ment and is therefore more easily combined with on-chip 
magnetic or electrostatic trapping techniques [27]. In- 
deed, based on these properties the storage of ensemble 
qubits encoded in collective spin excitations of an MDC 
has previously been proposed [28-30, 34 . However, the 
high densities and the much stronger coupling to phonons 
make the crystalline phase incompatible with conven- 
tional schemes for addressing and manipulating qubits 
encoded in individual molecules of the crystal. 

In this work we describe a new approach for QIP which 
is especially adapted to the physical conditions of an 
MDC. The main idea is to use the strong phonon cou- 
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pling to create a controlled phonon mediated interaction 
for gate operations between qubits instead of the direct 
dipole-dipole interaction. Similar gate operations that 
feature such a phonon bus are frequently employed in 
quantum computing and simulation schemes in chains of 
trapped ions [35H39] . Here light forces can be used to 
displace individual ions in a state selective way to cre- 
ate effective interactions between the otherwise highly 
decoupled internal states. We are interested in similar 
techniques by creating state dependent dipole moments 
for molecular qubits using microwave fields only. How- 
ever, as noted above opposed to the large spacings in ion 
chains the high densities of MDCs make it difficult to ad- 
dress single sites. To overcome this problem we suggest 
to combine the proposed technique with the concept of a 
marker qubit [40l|4T], i.e. an additional molecule which 
is trapped in a separated layer above the crystal. Due to 
strong inter-layer dipole-dipole interactions this marker 
qubit modifies the local phonon structure of the MDC 
such that phonon mediated qubit operations can be im- 
plemented at pre-specified lattice site in the crystal using 
otherwise global addressing only. As the marker moves 
between different designated sites of the MDC more ef- 
ficient gate operations between qubits located at distant 
sites of the crystal are possible in comparison to previous 
proposals relying on nearest neighbor interactions. 

The remainder of the paper is structured as follows. In 
Sec. |TI| we present a brief overview of MDCs and a qual- 
itative outline of the key ideas of this work. In Sec. Ill 
we describe in more detail the implementation of molec- 
ular qubits state with controllable dipole moments and 



discuss in Sec. [IV| the resulting direct and phonon me- 
diated dipole-dipole interactions. As a specific example 
we evaluate in Sec.|V|the implementation of a local gate 
operation in a ID MDC and finally summarize the main 
results and conclusions of this work in Sec. (VTl 




FIG. 1: A quantum register of molecular qubits confined in a 
2D plane and stabilized in a self assembled crystal by strong 
repulsive dipole-dipole interactions. The dipole moments of 
the molecules are aligned by an electric bias field and quan- 
tum information is encoded in long-lived spin or hyperfine 
levels of the rovibrational ground state. 



molecules is confined to the orthogonal plane by a strong 
transverse trapping potential (see Fig. [T]), dipole-dipole 
interaction between molecules are purely repulsive and 
at low temperatures the system is stabilized against close 
encounter collisions. The dynamics of such an ensemble 
of polar molecules is described by the Hamiltonian 
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where m is the mass and and the position and mo- 
mentum operators of the molecules and D = /i^/(47reo) 
for an induced dipole moment /i. At low temperatures 
and for strong transverse confinement the dynamics of 
the molecules is characterized by a dimensionless param- 
eter which is the ratio between interaction and kinetic 
energy for a mean interparticle spacing a, 



II. MOTIVATION AND OVERVIEW 

We describe the implementation of quantum gate op- 
erations in a ID or 2D MDC as shown in Fig. [l] where 
qubits are encoded in long-lived spin or hyperfine states 
of polar molecules and stabilized against collisions by 
strong repulsive dipole-dipole interactions. We start in 
this section with a brief review of the basic properties a 
MDC followed by a qualitative outline of the main ideas 
of our proposal. This overview then serves as a moti- 
vation and guideline for the more elaborate calculations 
detailed in the remaining sections of the paper. 



A. Molecular dipolar crystals 

Compared to neutral atoms the distinguishing feature 
of polar molecules is a large electric dipole moment of 
up to several Debye which can be aligned in the presence 
of an electric bias field. If in addition the motion of the 
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(2.2) 



At high densities the potential energy dominates and for 
r^i ^ 1 the system enters a crystalline phase j3TH33] , 
which for molecular dipole moments of several Debye oc- 
curs at typical interparticle spacings a of a few hundred 
nanometers up to ~ /im. 

In the crystalline regime molecules become localized 
at fixed equilibrium positions Tq, and, e.g, in 2D form a 
triangular lattice. The dynamics of the MDC is then well 
described by the residual small fiuctuations around the 
respective equilibrium positions, 2^ + Ui, which lead 
to collective phononic excitations of the crystal [30 and 



^MDC — hxJka\ak' 



(2.3) 



Here ujk and ak are the frequencies and bosonic annihi- 
lation operators for phonon mode k as we will discuss in 
more detail in Sec. IIVI 
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a) Rotational levels 

E 
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b) Fine and hyperfine 
substructure 




FIG. 2: a) Sketch of the rotational spectrum in the presence 
of an electric bias field and qubit states \g) and |e) are encoded 
in two long-lived rotational eigenstates of the molecule, b) 
Qubit states are encoded in two spin or hyperfine sublevels of 
the rotational ground state. Spin-rotation interactions ^ jsr 
allow the implementation of state dependent dipole-moments 
by admixing higher rotational states using microwave fields. 
See Sec. Illllfor a more detailed discussion. 



In contrast to trapped ions or neural atoms the state de- 
pendent dipole-dipole interactions naturally provide a di- 
rect qubit-qubit coupling between neighboring molecules. 
However, while the interaction (2.5) contains the basic 



ingredients for quantum information processing, it is in 
this form only of limited use for the implementation of 
quantum gate operations between individual qubits in 
a crystal. First, (2.5) describes global spin-spin interac- 



tions induced by homogenous static and microwave fields, 
and therefore would again require additional local manip- 
ulations on the small scale of the MDC lattice spacing 
a. Second, in writing Eq. (2.5) we have ignored fluctua- 



tions of the molecules around their equilibrium positions, 
which induce decoherence, if not appropriately taken into 
account. 



Resonantly enhanced phonon mediated 
dipole-dipole interactions 



B. Molecular qubits & spin-spin interactions 

The MDC represents a stable configuration where in 
reminiscence of a Wigner crystal of trapped ions [1 the 
molecules are well localized at individual sites of the self- 
assembled crystal lattice. We are interested in using this 
crystal as a quantum register where quantum information 
is encoded in two long-lived internal states \g) and |e) of 
the individual molecules. As illustrated in Fig. [2j rota- 
tional spectroscopy [42. allows us identify different sets 
of qubit states [22l[30], and for example \g) and |e) can 
be encoded in two rotational eigenstates or in two spin 
or hyperfine sub-levels of a single rotational manifold. 

While in principle both rotational and spin degrees 
of freedom would provide a suitable choice for long- 
lived qubit states we will focus below explicitly on spin 
qubits as shown in Fig. |2]b). In this case the idle qubit 
states have identical dipole moments while appropriate 
microwave fields can be used to admix higher rotational 
levels to engineer dipole moments in a state selective way. 



More precisely we will show in Sec. Ill how to generate a 
state dependent dipole moment which can be written as 



(2.4) 



Here (J.t 



-WzCFz is an arbitrary Pauli oper- 



ator in the qubit space {|^), |e)}, specified by the unit vec- 
tor w. To guarantee the stability of the crystal indepen- 
dently of the internal states we require e := d/j^/fio <C 1. 
Under these conditions the picture of molecules being 
localized around their equilibrium positions is 
still valid, but on top of that the state-dependent dipole- 
dipole coupling leads to additional direct spin-spin inter- 
actions of the form 
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(2.5) 



To overcome the limitations of direct spin-spin inter- 
actions mentioned above, we propose in this work the 
implementation of enhanced phonon mediated gate oper- 
ations which rely on virtual excitations of isolated vibra- 
tional modes in the crystal. To illustrate the basic idea 
of this approach we must now take into account mod- 
ulations of the dipole-dipole coupling induced by small 
position fluctuations of the molecules. As we show in 



more detail in Sec. |IV| these fluctuations result in spin- 
phonon interactions H^-p where for e <C 1 the dominant 
contribution is of the form 



i,k 
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(2.6) 



Here we have expressed displacements in terms of the 
normal mode operators such that ~ d/j. is the cou- 
pling of phonon mode k to the spin at lattice position 
i. 

To make use of the spin-phonon coupling in a con- 
trolled manner we introduce a modulation of the dipole 
moment Sfi{t) = Sficos{uJot) with a frequency uq. Un- 
der certain conditions detailed in Sec. IIVI we can then 
eliminate the external degrees of freedom and as a re- 
sult we obtain an additional, phonon mediated spin-spin 
coupling. 



Hi 




(2.7) 



We see that by tuning cuq close to a particular phonon 
frequency uJk we can in principle amplify the coupling 
mediated by that specific mode to the extend where it 
dominates over all other contributions, in particular also 
over direct spin-spin interactions in Eq. (2.5). In that 



case the spatial pattern of the resulting spin-spin inter- 
action is solely determined by the mode function of the 
resonantly enhanced phonon mode. 
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FIG. 3: Illustration of the basic steps of an entangling op- 
eration between the movable marker molecule and a register 
qubit in the crystal. See text for more details. 



D. Marker qubit s and local gate operations 

As a final ingredient for our proposal we now consider 
an additional "marker" qubit, which is independently 
trapped in a layer parallel to the crystal. By using a 
weak optical tweezer potential the marker molecule can 
be moved along the xy plane and locked to a specific reg- 
ister molecule in the MDC via strong attractive interlayer 
dipole-dipole forces. This results in the desired forma- 
tion of a localized phonon mode which ideally involves 
only a single molecule in the crystal. Thus, by employ- 
ing the resonantly enhanced, phonon-mediated spin-spin 
coupling technique described above, this allows us to im- 
plement a general set local two qubit operation between 
the marker qubit and a target register molecule r, 

ZY£(tZl,(/.)=e^^----. (2.8) 

Note that this set of operations is universal and includes 
single qubit rotations of the register molecule by applying 
UQ{w^(j)) with the marker qubit prepared in one of the 
eigenstates of a^. 

As illustrated in Fig. |3]the implementation of arbitrary 
quantum information processing protocols in a MDC can 
then be decomposed into the following steps: 

1. The trapping layers for the MDC and the marker 
molecule are well separated such that the state 



of the marker qubit can be manipulated indepen- 
dently of the register qubits in the crystal. 

2. By using a weak optical tweezer potential Vtw(^, y) 
the marker qubit is positioned in the xy plane above 
a target memory qubit r and the optical barrier in 
z-direction is reduced adiabatically. 

3. The marker and the register molecule r form a lo- 
calized phonon mode which is used to implement 
a local phonon-mediated gate operation UQ{w^(j)) 
while leaving the rest of the quantum register un- 
affected. 

4. The two planes are separated again and the marker 
molecule can be guided to another lattice site s 
where the sequence is repeated to generate entan- 
glement between the two memory qubits. 

Different generalizations of this basics scheme using mul- 
tiple marker molecules or multiple layer configuration can 
be envisioned and would allow for a speed up and a par- 
allelization of QIP protocols. However, the key ingredi- 
ent in all those schemes is the implementation of local 
o'wd-uj operations and in the remainder of this work we 
now present a more detailed discussion of how this can 
be achieved in the context of polar molecules. 

III. MOLECULAR QUBITS AND STATE 
DEPENDENT DIPOLE MOMENTS 

In this section we start with a discussion of molecular 
qubit states and show how we can use combinations of 
static and microwave fields to generate the desired dipole- 
dipole interactions between two molecules. As the inter- 
nal structure of molecules can vary significantly between 
different species we will here illustrate the implementa- 
tion of our scheme for the case of a ^S-molecule, but 
also outline the general strategy for achieving a similar 
controllability for other molecules. 

A. Molecular spin qubits 

With the electronic and vibrational degrees of freedom 
frozen out, the internal dynamics of cold molecules is 
well described by its rotational spectrum. For S ground 
state molecules this spectrum is usually dominated by 
the rotational energy, which in the presence of an elec- 
tric bias field E5 = Ei^e^ is given by i^^ot = — (izEi). 
Here N is the angular momentum of the nuclei, B is 
the rotational constant which is typically of the order 
of 10 GHz and fiz = /J^b'f^z is the z-component of the 
molecular dipole operator. The bare dipole moment jUb 
depends on the molecular species and can be as large 
as 8.9 Debye for SrO. In the limit ^5^0 the eigen- 
states i^rot are the usual angular momentum eigenstates 
\N^Mn)o and we obtain the anharmonic energy spec- 
trum En^Mn — B{N + 1)N of a rigid rotor. At finite. 
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a) Rotational Spectrum 



b) Dipole moments 
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FIG. 4: Rotational spectrum in the presence of an electric 
bias field E^. a) Energy levels of ^rot = BN^ — jlzEh as 
a function of the bias field Eh^ where curves for eigenstates 
\N^Mn)ei, are labeled by (A^, |MAr|). b) The corresponding 
induced dipole moment /i = Ei,{N, MN\fiz\N, Mn)ei, in units 
of the bare dipole moment of the molecule /j^b- 



but not too large this level structure is still approx- 
imately correct, but the new eigenstates \N^MN)Eb 
now superpositions of states with different N values and 
therefore exhibit a non- vanishing induced dipole moment 
eA^^ ^N\i^z\N, Mn) Et, 7^ 0. Fig. [4]shows the energies 
and induced dipole moments for the lowest rotational 
states as a function of the applied electric bias field. We 
see that in principle we can select two of these levels as 
our qubit state \g) and |e) and by adjusting the bias field 
induce similar or state dependent dipole moments. How- 
ever, this would only provide a limited amount of control 
and for more flexibility we focus in the following on qubits 
encoded in spin sublevels of a single rotational state. 

For molecules with an unpaired electron the rotational 
spectrum is in a good approximation described by the 
Hamiltonian 1^ 



i^rot = - jl.Eb + 7srNS + i^hyp, 



(3.1) 



where S is the electronic spin and 7sr the spin-rotation 
coupling constant. The Hamiltonian ^^hyp accounts for 
additional hyperfine interactions, which in the presence 
of a non-zero nuclear spin / 7^ lead to a more involved 
level structure as described below. However, this addi- 
tional spin structure does not affect the main arguments 
presented here and for simplicity we can assume / = 
for the moment. The coupling between spin and rota- 
tion ^100 MHz is typically much weaker then the 
rotational energies and for a sufficiently large bias field 
El) ~ B / jii rotational and spin degrees of freedom are to a 
good app roximation decoupled. Therefore, eigenstates of 
Eq. ([3^ can be written as ^ \N,Mn)e, \S,Ms), 
and, e.g., for S = 1/2 each rotational level will simply 
split into a spin doublet. In particular, we now choose 
two spin sublevels of the rotational ground state for our 
qubit states l^f) and |e) as shown in Fig. |2j In this case 
both states have identical induced dipole moments, which 
makes these state highly insensitive to noisy electric fields 
[3Qj . Nevertheless, as we will describe in the following, we 
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FIG. 5: General scheme for the creation of a state dependent 
dipole moment by admixing a rotationally excited level |r) to 
the qubit states \g) and |e). 



can use microwave fields to couple \g) and |e) to higher 
rotational states where a finite admixing of order %r/B 
between spin and rotation allows us to induce state de- 
pend dipole moments. 



B. State dependent dipole moments 

We now outline a general strategy t o gen erate state de- 
pendent dipole moments of the form (2.4). Restricted to 



the qubit subspace {|^), |e)} an arbitrary dipole operator 
is given by 



(3.2) 



where /i^^ = {u\fi\v) and \u), \v) e {\g), |e)}. By encod- 
ing qubit states in two spin sublevels of the rotational 
ground state, we obtain /j,gg = 11^^ = /j.^ a nd /j^^g = and 
the internal state dependent part of (3.2) disappears for 



such an idle qubit. To introduce non-trivial interactions 
the states \g) and |e) are coupled by microwave fields [44] 
to an auxiliary, rotationally excited state |r) as shown in 
Fig. [sj similar to [39 . This auxiliary level is character- 
ized by an induced dipole moment /j,^^ ^ fi^ and small 
transition dipole matrix elements fi^g = fi^^ <C /[Xo^Mrr 
to eventually suppress unwanted flip-flop interactions as 
discussed below. We will show in detail how such a three 
level system can be realized in the following sections. 

In the lambda configuration shown in Fig. |5] the states 
1^) and |e) are coupled to a single rotationally excited 
state |r) with respective Rabi frequencies ftgc'^^a and 
l^ee*"^^ and a detuning A. We introduce the ratio of the 
Rabi frequencies tan(i^) = Qg/fte and (^eff = \/^9 ~^ ^e- 
The difference between the phases is denoted by 2r] = 
(j)g — (j)e. In the frame rotating with the microwave fre- 
quency cjmw = cj^ — a, the three eigenstates of this sys- 
tem are given by 

\a^) = cos(^c.)|ao) + sin(^c.)k), (3.3) 
= |/3o), (3.4) 
Ijq) = cos(^^)|ao) + sin(^^)|r), (3.5) 
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where 2tan(^c.) = -(A/lleff) + V'4 + (A/l^eff)^ and 
2tan(^^) = (A/l]eff) + \/4 + (A/l^eff)^- Astan(^c.) tends 
to zero when the effective interaction strength disap- 
pears, the states and \/3q) are adiabaticahy con- 
nected to the superposition states 



\ao) = sin(i^)e*^|^) + cos(i^)e 



(3.6) 
(3.7) 



which correspond to the so-cahed bright and dark states 
of a given set of microwave parameters. 

The original qubit states can be expressed in terms of 
the dressed eigenstates, 

1^) = e-^^[sin(^)|ao)+cos(^)|/3o)], (3.8) 
|e) = e^^[cos(^)|ao)-sin(^)|/3o)], (3.9) 



which shows how our qubit states are modified when the 
fields Clg and are adiabatically turned on. We can now 
determine the modified dipole matrix elements in (3.2) in 
the presence of the fields as 

A = Mol + 2AXi|ao)(ao|, (3.10) 

where we have defined the dipole moment 

Ml = ^ sin^(^«)(Mrr - Mo)- (3.11) 

By introducing a normalized vector w = 
(sin(2i^) cos(27^), sin(2i^) sin(27^), cos(2i^)) the projec- 
tor |ao)(ci^o| can be expressed in terms of the usual Pauli 
matrices for the qubit states \g) and |e) as 



2|ao)(ao| = 1 



(3.12) 



The complete control over our system parameters thus 
allows us to create an arbitrary state dependent moment 
IJiiCFy^ by adiabatically turning on the microwave fields. 
The important timescale at which the AC-fields can be 
varied must be much longer than 1/A, i.e the time span 
r it takes to turn on the microwave fields to strength l^eff 
must fulfin r > l^eff/A^ 



Finally, we remark that result (3.10) differs slightly 
from the purely state dependent moment promised in 
(2.4) as we generate the projector |ao)(c»^o| rather than 
cFy^. Also we obtain either a positive or a negative dipole 
moment /l^^, while ideally we would be interested in gen- 
erating oscillating dipole moments /L^i(t) ~ cos(co'ot) with 
zero mean. These fiaws can easily be corrected by using 
first of all a time depend bias field ^h(f) such that the 
state independent part is absorbed, AXo(^) + Mi(^) = Mo- 
This requires a similar adiabatic condition as above. Sec- 
ond, to obtain purely oscillating fields we point out that 
via the relation = — we can change the sign of the 
state dependent dipole moment by simply changing the 
microwave field control parameters v and 77. Therefore, 
within the adiabatic regime defined above, the three level 



configuration shown in Fig. [5] is sufficient to engineer ar- 
bitrary time dependent dipole moments of the form 



jl = /Lfcol +Ml(^)^^ 



(3.13) 



The resulting dipole-dipole interaction between two 
molecules then becomes 



• MoMi(^) + ^ii) + Mi(^)'^L-^^, (3.14) 



which we will use as a starting point for our discussions 
of effective spin-spin interacti ons i n Sec. [IV| below. 
Note that the form of Eq. (3.10) and in particular the 



validity of Eq. (3.14) are based on the important condi- 
tion /j^rg = Mre ^ 0. If this is uot the case additional cor- 
rections to V^^ appear, for example in form of a ffip-ffop 



interaction scaling as tan^(^cK)|Mer P + To 

avoid such unwanted contributions we must identify a set 
of states 1^), |e) and |r) where dipole transition matrix 
elements are strongly suppressed. In the following we de- 
scribe in more detail in two examples, how this can be 
achieved in practice. 



C. Example I: electronic spin qubits 

We first consider the case of a molecules with 
5* = 1/2 and zero nuclear spin. The lowest states of the 
resulting rotational spectrum are shown in Fig. [6] for a 
typical value of the spin rotation coupling %r/B ^ 0.01. 
In the following we label eigentstates by \N^Mn^Ms) 
with the index £^5 omitted and we encode our qubit 
in the two spin states of the rotational ground state, 
1^) = |0, 0,-1/2) and |e) = |0, 0,1/2). In the idle case 
these two states have identical induced dipole moments, a 
zero transition dipole moment and are well separated by 
the rotational splitting ~ 2B from the higher rotational 
states. 

We identify the state |2,2, -1/2) in the A/" = 2 mani- 
fold as a suitable level |r) that has zero or spin forbidden 
transition matrix elements with both \g) and |e), while 
its induced dipole moment fij-r differs significantly from 
the one in the ground state. The state |e) can be coupled 
directly to |r) with a high ac-field amplitude to overcome 
the strongly suppressed transition dipole matrix element. 
For the level \g) we must use a two photon process via 
the intermediate level \x) = |1, 1,-1/2). The proposed 
microwave fields also couple to additional unwanted lev- 
els 1^) = |1, 1,1/2), \s) = |2,2,l/2) and |zi,2) of our 
spectrum, as outlined in Fig. |6] 

We are interested in a regime where all contributions 
from levels other than |^), |e) and |r) are small. Here 
the spin-rotation coupling plays a key role. It allows 
transitions between opposite spin states while at the same 
time it suppresses unwanted transition matrix elements 
of spin forbidden transitions by Jsr/B. In addition 
lifts the degeneracy between states \x) and \y) as well as 
|r) and \s) and enables the design of a state dependent 
coupling to higher levels. 
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FIG. 6: Schematic plot of the level structure of molecules. 
We make use of the specific features of this model resulting 
from small spin-rotational couplings "jsr/B <^ 1 to create an 
effective A-system between the levels \g), \e) and |r) by ap- 
plying the microwave fields Oi, Q2 and Oe- 



For 5 ^ r^i,r^2,A we can adiabatically eliminate \x) 
in the two photon process to give an effective Rabi fre- 
quency Qg = ^^1^2/^ where remaining contributions of 
the unwanted level are estimated to be of the order of 
Qi/5 and ^2/^ or smaller. The transition dipole elements 
between |e), \y) and \s) are similar to the ones between 
1^), \x) and |r). By making use of the splitting between 
between the opposite spin states \x) and \y) as well as |r) 
and \s) the two unwanted levels can be eliminated if the 
spin rotational coupling fulfills 7sr > S. Then the leading 
contributions from \y) and \s) are of the order of Qi/S 
and QiQ2/^'^ respectively. The ac field that couples |e) 
directly to the opposite spin state |r) induces additional 
spin allowed interactions between \g/e) and 1 2^1/2) with 
large Rabi frequencies of the order of /jrs)- How- 
ever the detuning to these levels, ujz — oJs is of the order 
of the rotational constant such that these levels admix 
only in orders of l^g/Trs- 

Let us now recollect all conditions and parameters 
involved. We require a small spin rotational coupling 
Jsr/B ^ 0.01 to strongly suppress transition matrix el- 
ements between |e) and |r). To eliminate the levels \y) 
and \s) the detuning S must not exceed the spin rotation 
coupling, S ^ and the Rabi frequencies must fulfill 
^1,2 S. The detuning A should be of the same or- 
der as Qg^e to obtain reasonable admixings of the upper 
level. They are then put into relation with other energies 
by Qg ~ QiQ2/^- As the rotational constant is typically 
of the order of a few to tens of GHz and 7sr ~ 100 MHz we 
obtain an upper bound for rti^2 around 10 MHz and effec- 
tive parameters Qg^^^^ ^ 1 MHz. Therefore, the lowest 
internal energy splitting A in this model still exceeds the 



relevant dipole-dipole interaction energy Vdd < 100 kHz 
such that the influence of other molecules can be ne- 
glected. 



D. Example II: nuclear spin qubits 

As a second example we now consid er m olecules with 
a non-zero nuclear spin, which in Eq. (3.1) leads to ad- 



ditional hyperfine interactions of the form 
^hyp = ^SI + c(Sn)(ni). 



(3.15) 



Here n is the unit vector along the internuclear axis and 
b and c are the hyperfine coupling constant which typ- 
ically are much smaller than but comparable to the 
spin-rotation coupling, 7sr ~ l^hypl ~ 27r x 100 MHz. 
For the simplest case 5 = / = 1/2 we obtain the level 
diagram shown in Fig. ^ In particular, we see that the 
rotational ground state splits into four hyperfine levels 
corresponding to the eigenstates of the total spin oper- 
ator J = S -\- I and the same structure as above re- 
peats itself in the = I^Mn = 0) rotational manifold. 
Therefore, we can choose the states |^) = |7V = 0, Mn = 
0, J = 1, Mj = —1) and |e) = |0, 0, 1, 1) as our qubit and 
the auxiliary, rotationally excited state |r) = |1, 0,0,0). 
In contrast to the spin qubit discussed above, both \g) 
and |e) are now coupled to |r) via spin- forbidden tran- 
sitions with strongly suppressed dipole matrix elements 

/ier — l-f^gr 0. 

Similarly to the previous model there are various spin 
allowed transitions between the = an the = 1 
manifold that lead to resonant flip-flop processes of the 
form |r)i|e)2 \y)i\x)2 and drive the system outside of 
our three level subspace. These processes can be ener- 
getically suppressed by introducing additional ac-Stark 
shifts as indicated in Fig. [7| 



IV. PHONON MEDIATED SPIN-SPIN 
INTERACTIONS IN A DIPOLAR CRYSTAL 

In the previous section we have shown how microwave 
fields can be used to manipulate the dipole-dipole inter- 
action between two molecules. We now extend this dis- 
cussion to molecules embedded in a MDC where state de- 
pendent dipolar forces introduce additional interactions 
between the external and internal degrees of freedom. 
Our main goal is then to show how time dependent con- 
trol fields can resonantly enhance the coupling to a par- 
ticular phonon mode and lead to controllable effective 
qubit-qubit interactions after elimination of the motional 
degrees of freedom. 

The complete Hamiltonian of our model is given by 



tot 



E 



Pi 

2m 



^;ot + Hrap(r.) +Vdd({rJ), (4.1) 
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Mn + Mj 



FIG. 7: Schematic level structure of molecules with a nu- 
clear spin / = 1/2. The qubit states \g) and |e) are coupled 
to the auxiliary state |r) via microwave fields of Rabi fre- 
quencies Qg and Qe- As described in the text, the microwave 
field Qs is used to introduce additional ac-Stark shifts for the 
levels \x) and \x) to suppress otherwise energetically allowed 
transitions into those states. 



is in general given by 

^2 



MDC 



D ^ 1 -3cos(i9i^) 



••|3 



(4.5) 

where 6ij is the angle between — and the z axis. The 
second term contains internal degrees of freedom only 
and is given by 



4 = 5]4'ot+V^dd({r?}). 



(4.6) 



Note that a residual spin dynamics from H^^^ can be ab- 
sorbed by applying additional global 5-fields and we can 
neglect this term in the following discussion. Finally, the 
last term in Eq. (4.4) contains all remaining contribu- 



tions to the dipole-dipole interaction which affect both 
internal and external degrees of freedom. It is defined 
as i^s-p = Md({^«}) ~ Md({^i }) leads to non-trivial 
spin-phonon interactions, which will be the main focus of 
our discussion below. 



which includes the kinetic and internal energy molecules, 
the external trapping potential and the full dipole-dipole 
interaction 

,y .r 1 ^ f^i-P^j -^i^ij-Mi^ij-P^i) 



87reo ^ 



.,13 



(4.2) 

with Uij rij/\rij\ and Vij = - Vj. 

Under conditions detailed in the previous section we 
can restrict our discussion to dipole operators of the form 



(4.3) 



with all dipole moments aligned along the z axis. It is 
then convenient to split the dipole-dipole interaction into 
two parts as VddiUi}) = "KdKi^^J) + ^ddiUi})^ where 
V^{{Ti}) (X /iQ depends only on the molecular positions 
while V^dii^i}) ^ S/i^S/j? contains the internal state de- 
pendent part of the interaction. We are interested in 
the regime where a stable crystalline phase is realized 
independently of the internal configuration which is the 
case when the stabilizing force V^{{ri}) dominates over 
Vj^{{Ti}). Then the molecular position operators are 
given by 



A. Phonons 

The basis of our description is the formation of a molec- 
ular crystal in our model described by the crystal Hamil- 
tonian 

^MDC = ^ + Hrap(ri) + ^ ^dd(ri - r^), (4.7) 



with Vdd(i*) = D{1 - 3nl)/r^ and D = /Hq/Attco. In the 
crystalline phase we can follow the standard procedure 
j43] and expand the potential to second order in the small 
displacements around equilibrium positions that are 
determined by minimizing the potential energy in Eq. 
(4.7). While in homogeneous ID or 2D systems the 



would form a regular lattice in the general case details 
depend on the external trapping potentials as described 
m more detail in Sec. [V| The expansion of Hamiltonian 
(4.7) is given by 



MDC 



^ 2m 



9rr, 9 



(4.8) 



where = 5,«a^/3 [Fdd({r?}) + Hrap({r?})]. To diag- 

f'.^'"" ^\ -''^—''- '^--^- --^--^ onalize the quadratic Hamiltonian we introduce a set of 

tions around the equilibrmm positions and the total . 



u,- where u,- denote small fluctua- 



Hamiltonian can be decomposed into three parts. 



bosonic mode operators 
placements Ui^a by 



that are related to the dis- 



MDC 



(4.4) 



Here the first term represents the external molecular mo- 
tion and describes the formation of the crystal and its 
vibrational excitations. This term has already been in- 



(4.9) 



troduced in Eq. (2.1) for a purely 2D configuration and 



The normalized modefunctions ^ind eigenfrequencies 
ujk are obtained from the resulting eigenvalue problem 



9 



and we can write the crystal Hamiltonian in terms of the 
bosonic operators as 



C. Effective spin-spin interactions using oscillating 
dipole moments 



^MDC ^ fujJka\ak' 

k 



(4.10) 



B. Spin-spin and spin-phonon interactions 

In presence of the control fields the dipole-dipole in- 
teraction V^/d({^*}) giv^s additional molecular couplings 
that strongly depend on the internal degrees of freedom. 
To maintain the stable crystalline phase that has been 
described previously such additional interactions must 
be small compared to the stabilizing potential V^^{{Yi}) 
such that the molecular displacements from the original 
equilibrium positions fulfill |(ui)| <C a for arbitrary in- 
ternal configurations. This is realized when 



Mo 



e < 1. 



(4.11) 



The state dependent part of the dipole-dipole interaction 
is given by 



Up to this point ^g-p describes an uncontrolled cou- 
pling of the spins to the crystal phonons. It is our aim 
to amplify the interaction with modes around a chosen 
frequency ujq. To achieve this we let the small dipole 
moment dfi oscillate as 



5ii{t) = cos(co'ot), 



(4.17) 



with the picture of a classical harmonic oscillator that is 
driven close to its resonance frequency in mind. This can 
be achieved by time-varying control parameters under 



conditions detailed in Sec. |IIIB[ We proceed in a rotating 
frame with respect to H^^^ = ^^hujQa\ak^ where the 
total Hamiltonian (4.4) takes the form 



^tot =^h/^ka\ak + ecos(a;ot)^Vdd(r?^)aJ^ 

k iy^j 



1 + cos(2a;o^) 



(4.12) 



and can be expanded to leading order in u^. The lowest 
order of the expansion gives the spin Hamiltonian (4.6) 



where we have introduced an effective magnetic field 
^eff = ^j^i ~ r^p. The remaining part of the 

expansion leads to spin-phonon interactions, and to first 
oder in u.- we find 



6 



k,i 



Here we have introduced 



(4.14) 



(4.15) 



where = ^fe derivative of the dipole- 

dipole potential is 



3^ / n^(5n2 - 1) 
^ddW = ^yi^^l - 1) 



(4.16) 



Here the A/^ = cj^ — co'o denote the detunings between 
ujQ and the phonon frequencies ujk. Note that we have 
neglected terms of the order O(e^u^) that arise from the 
second term on the righthand side of Eq. (4.14). Eventu- 
ally such contributions only give corrections of the order 
of O(e^), and - as they only couple to low frequency 
modes - they are not enhanced by the resonance condi- 
tions discussed in the following. By making a rotating 
wave approximation we obtain 



^tot - ^ hAkCilak 



^^^dd(r?,)^j^^^ 



^Atai(al + afe), (4.18) 



kA 



which is valid for 



^dd 



(4.19) 



This essentially limits our choice of cuq to frequencies that 
are of the order of the dipole-dipole interaction D/a^ or 
above. 



Based on the effective Hamiltonian (4.18) we now elim- 



inate the phonons by a canonical transformation. This 
Lang-Firsov or polaron transformation [421 145] is defined 
by H ^ (jHtj^ where 



tj = e^ with S = Y^ 



-aj,(a,-4). (4.20) 



This transformation changes into a picture where internal 
and external degrees of freedom are decoupled as Htot = 



10 



MDC 



s,eff 



Ep , where 



H. 



s,eff 



^2 . ^2 



(4.21) 

The phonon couphng has thus been ehminated in favor of 
an additional spin-spin interaction, the phonon mediated 
interaction (PMI) Vpm, and a constant energyshift Ep = 
^ Xl'^e^ /2hAk known as the polaron shift. The first 



term in (4.21) is what remains of the direct spin-spin 



interaction. 



D. Discussion 

The phonon mediated interaction can be enhanced by 
reducing the detuning for a certain phonon mode. 
This is of course only valid up to a certain point, as we 
can see from the additional collective molecular displace- 
ments 



UatU^ =at 



^ 2hAi 



(4.22) 



We introduce an upper bound to the relative displace- 
ment for any internal configuration as Au = max{|u^ — 
u^+i |}. Then the phonon mediated interaction can be en- 
hanced as long as Au does not exceed a certain fraction 
of the mean interparticle spacing. 



A" = I E(cl - cl 



+1^ 



2muJk 



E 



hAk 



<Ca. (4.23) 



This condition connects e and A/, such that smaller dipole 
ratios e allow us to choose cjq closer to resonance with a 
particular phonon mode of frequency uor^ as illustrated in 
Fig. [sj)) for a toy model specified below. For a spectrum 
that is sufficiently sparse around the resonant mode, we 
can achieve A^^ <C A/e, \/ k ^ R such that the phonon 
mediated interaction is to a good approximation given by 



pm 



2 ^ 



hAi 



(4.24) 



As the sign of A^^ can be freely chosen by tuning cjq above 
or below the respective resonance frequency uor the PMI 
can be made attractive or repulsive. By selecting a fixed 
value for Au we connect e and A^^ such that the PMI 
becomes effectively of first order in e. In this manner the 
ratio of the direct interaction over the phonon mediated 
one can be made arbitrarily small by variation of the 
dipole ratio which of course comes at the cost of reducing 
the absolute value of the interaction strength. 

Other properties of the PMI are revealed when study- 
ing A*^ (X Ej^dd(i*ij)(Cfc - Ci) which is nonzero only 
if there is relative motion between molecule i and any 
other molecule j in the direction of a nonzero gradient 



a) Two molecules in a harmonic trap 




b) Phonon mediated interaction 
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FIG. 8: a) Phonon mediated interactions between two polar 
molecules confined to the x-axis by strong perpendicular con- 
finement and with dipole moments aligned along the z axis. 
The system is stabilized in x-direction by a harmonic trap- 
ping potential v ^ v^. b) The strength Vpm of the phonon 
mediated interaction in this system with cjq tuned close to the 
frequency of the longitudinal breathing mode ujr — \f^v. In 
red (green) we outline the forbidden parameter regime for cjq 
when the molecular displacements are restricted to Au[e) <^ 1 
for different dipole ratios e. Note that the direct spin-spin in- 
teraction in this model has the value of 1/2 in units of e^D/a^. 



of Vci(i{r^j) when oscillating in mode k. Therefore, dom- 
inant contributions to the phonon mediated interaction 
arise from the high (low) energy longitudinal (perpendic- 
ular) modes in MDC's. In addition, every term in the 
sum is suppressed by a factor l/(r^^)^, which generally 
leads to a fast decay of PMI. 

In the following we brieffy describe a simple system of 
two molecules in a harmonic trap as shown in Fig. |8^) , 
which illustrates the basics features of our model. The 
molecules with dipole moments aligned along the z-axis 
are conffned to the x-axis by a strong transversal trap- 
ping and Vj.*^ap ~ mu'j^{yf-\-zf)/2-\-mu'^Xi/2^ with u± > u. 
The equilibrium spacing along the x-axis is given by 
a = (6/io/^^^)^^^- III this setup the oscillations in x, y 
and z directions decouple with two modes in every direc- 
tion. They include a center of mass (COM) mode at the 
respective trapping frequencies and in each direction a 
breathing mode with respective frequencies Ux/i' = >/5, 
^y/^± = \/ 1 — i^/i^± and uJz/^± = \/l — Zv jv^. The 
eigenvectors in every direction are simply (1,1)/V^ for 
the COM mode and (1, — 1)/>/2 for the breathing mode. 

The COM modes do not contribute to the PMI as 
there is no relative motion between the molecules. Also 
the perpendicular breathing modes do not contribute be- 
cause in our setup the ffeld gradient in z and y directions 
disappears along the x-axis. The only mode that medi- 
ates the PMI is the longitudinal breathing mode and the 
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MDC 



FIG. 9: Schematic setup for a double layer configuration 
where an optical standing wave potential confines the MDC 
and a single marker molecule in two planes separated in z- 
direction by h — A/2. The MDC has a lattice spacing a and 
for dipole moments aligned along generates an attractive 
periodic potential Vmdc for the marker molecule. For 6/a < 1 
the marker molecule locks onto the target register molecule 
and modifies the local phonon structure. 



resulting PMI strength is plotted in Fig. [sj)). By opti- 
mizing the detuning for a fixed value of ISaija = AiZ we 
obtain 



enh 



(4.25) 



As pointed out above this enhanced interaction depends 
linearly on the ratio e in contrast to the scaling of 
direct spin-spin interactions. 



V. MARKER MOLECULES AND LOCAL GATE 
OPERATIONS IN A ID CRYSTAL 

In the previous section we have described how we can 
resonantly enhance phonon-mediated dipole-dipole inter- 
actions by tuning the oscillation frequency cjq close to the 
frequency of a particular mode phonon mode. We now 
want to discuss in more detail how this technique can 
actually be exploited to generate local gate operations 
between qubits in a large dipolar crystal. To do so we 
focus on the setting shown in Fig. [9j In this double layer 
configuration a single "marker" molecule located in the 
upper trapping layer can be moved in the xy plane by 
a weak optical tweezer potential V^^(x^y). By bringing 
both layers close together, the marker molecule will lock 
to one of the molecules in the crystal and thereby create a 
spatially localized phonon mode. We can use this mode 
to implement a local gate between the marker and the 
target register molecule. Note that while below we will 
for simplicity only discuss a ID model, all arguments and 
conclusions presented in this section are equally valid for 
the 2D case. 



A. Local phonon modes of the bilayer system 

We first consider the case where the two optical trap- 
ping layers are brought close together and the tweezer 



potential is switched off. The minimal layer separa- 
tion h is then fixed by the wavelength of the trapping 
laser, h = A/2. Within the first layer polar molecules 
with dipole moments aligned perpendicular to the trap- 
ping plane experience repulsive dipole-dipole interactions 
and form an MDC with average interparticle distance a. 
The marker molecule resides in the second layer and sees 
a periodic potential Vmbc generated by the crystalline 
structure. Since this interaction is attractive the marker 
molecule 'locks' onto one of the molecules in the lattice 
below. This configuration has been analyzed in detail 
in Ref. [36] where it has been shown that the tunneling 
rate of such particles to neighboring sites is exponentially 
suppressed for rd^ I and b/a < 1 [46] . 

The presence of the additional marker molecule dis- 
torts the lattice and modifies the phonon spectrum. How- 
ever, already simple estimates show that down to values 
of 6/a ~ 0.5 the influence of the marker molecules on the 
lattice structure as well as the crystal phonon spectrum is 
relatively small due to the anisotropy and the fast decay 
of the dipole-dipole coupling. Therefore, as a main effect 
we expect the emergence of localized phonon modes due 
to the strong coupling between the marker and a single 
register molecule. 

As an example we plot in Fig. 



10 a) the phonon spec- 



trum of a ID dipolar crystal in the presence of the marker 
molecule. For this plot we have evaluated the modified 
equilibrium positions for this trapping configuration and 
solved numerically for the eigenvalues ujk for a system of 
= 50 particles assuming periodic boundary conditions. 
The acoustic phonon branch corresponds to longitudi- 
nal vibrations of the lattice and the two optical phonon 
branches describe transversal excitations with frequen- 
cies that are mainly determined by the transverse con- 
finement potential, i.e., cj/e^^ ~ uj±. We further observe 
the appearance of three modes ^^x^y^z with frequencies 
that are separated from the rest of the spectrum and 
correspond to isolated phonon modes which are localized 
around the marker molecule. This is illustrated more 



clearly in Fig. 10 3) where we plot the components of the 
modefunction Q. for the modes ^^x^y^z- we bring the 
layers closer together these modes become more and more 
localized and eventually reduce to relative oscillations be- 
tween the marker and the target register molecule only. 



B. Implementing local gate operations 

The ability to generate a spatially localized and en- 
ergetically separated phonon mode allows us to use 
the phonon mediated coupling technique described in 
Sec. |IVD| to implement a controlled gate operation be- 
tween the marker and the target register molecule. To do 
so the modulation frequency cjq is tuned close to one of 
the local mode frequencies, e.g uj^^^ to enhance the local 
PMI over the global direct spin-spin coupling. Since in 
general both interactions will be present, to discuss the 
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FIG. 10: a) Phonon spectrum for a homogenous ID MDC in 
presence of a marker molecule as shown in Fig.pfor a perpen- 
dicular trapping frequency uj± = ^^D/a^ /r]/ . Oscillations 
in different directions decouple and result in two optical and 
one acoustical branch. This is shown in the presence of the 
marker (blue, dashed) for h/a — 0.8 and without the marker 



(gray, solid). The three local modes uj^^ 
displayed in orange (red) for h/a — 0.9 (O.i 



and cji°" are 
b) The localiza- 



tion of ( 



is visualized by plotting mode functions of the 



marker (purple), the target molecule (green) and the max- 
imum of all others (brown). The dotted, solid and dashed 
lines represent values for the different modes cj^^^, u}y^ and 
cjj^^ respectively and the red area outlines the limit of the 
perpendicular stability for this setup 46 . 



fidelity of the local gate operations we write the effective 
spin Hamiltonian (4.21) as 



4,efr = t/o^S^L- + ^ E ' ^^.^j^^i-. (5.1) 

where the sum runs over all molecule pairs, excluding the 
the index pairs {m^r). Our goal is to identify the condi- 
tions under which Uq ^ Ures where Ures •= max{|/7ij|} is 
the strength of unwanted residual spin-spin interactions. 

We now focus on the setting described above where the 
marker molecule locks onto the register molecule such 
that rirnr = ^z- Then the only mode that contributes 



to the local interaction is 



,loc 



as oscillations in different 



directions decouple in this configuration similarly to the 
example given in Sec. |IVD| This mode has the additional 
advantage that up to a first approximation it does not 
contribute to the interaction between other molecules in 
the crystal. 

We can derive analytic approximations for the re- 
sulting interactions under the assumptions that i) only 
marker and pinned molecule oscillate with all other 
molecules at rest and ii) that the marker does not signif- 
icantly perturb the crystal lattice with molecular equi- 
librium positions exactly at the center of the transverse 
trap. This allows us to write the total local interaction 
approximately as 



Uo 



D 



±3- 



a3 \{b/a)^ {h/a) 



(5.2) 



a) Model Interactions b) Local VS Global Interactions 
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FIG. 11: a) Different contributions of direct and enhanced 
phonon mediated intera ctio ns to [/q and [/res for a similar sys- 



tem as outlined in Fig. 10 



The PMI is enhanced about cjf 
and gray dotted lines show the respective analytic approxi- 
mations. A contour plot of the ratio Uq /lives is given in b) for 
the specific choice of maximal displacement Au/a = 0.1. For 
comparison the absolute value of Uq is plotted in the same 
figure to illustrate the regime of interest. 



Here the two contributions arise form the direct coupling 
and the enhanced phonon-mediated interaction respec- 
tively for a specific choice of Au = Au/a respectively. 
Under the similar assumptions we can show that the 
magnitude of the residual global spin-spin interactions 
is given by 



3 2{b/af -S{b/a) 



4 (1 + (5/a)2)7/2 



(5.3) 



which is dominated by {D/a^)e^ /2 for reasonable values 
of b/a and e. Exact numerical calculations for the differ- 
ent contributions to the local and global interactions are 
summarized in Fig. 11a) and show good agreement with 



our analytic approximations. 

As one of the main results of this work Fig.[lT]b) shows 
the relative magnitudes of local and global spin-spin cou- 
plings for different system parameters e and b/a. For this 
plot we have enhanced the phonon mediated interaction 
up to the limit Au/a = 0.1 and have chosen similar signs 
for direct and phonon mediated contributions. Under 
these conditions we see that the local coupling can by far 
exceed residual direct spin-spin interactions. 

In summary we have shown that by employing local 
modes in our effective spin-spin model the local couplings 
can be amplified over the residual global interactions in 
the crystal. While in principle the ratio Uo/Ures ~ 1/e 
can be made arbitrary large it is in practice limited by 
the increase of the overall gate time. Nevertheless, ra- 
tions of Uo/Ures ^ 100 can already be obtained for rea- 
sonable values of e < 0.1 which still allows us to exploit 
more than 1/lOth of the dipole-dipole interaction D/a^ 
for gate operations. The local PMI is enhanced by the 
local direct coupling ~ l/{b/a)^ but even for large dipole 
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ratio s e ^ Au it exceeds the latter, as can be seen from 
( |5.2[ ), and shows the advantage of our model over purely 
direct dipole-dipole interaction. 



C. Implement ion using optically trapped molecules 

While so far our discussion of the gate operation has 
been quite general we finally want to address several is- 
sues specific for an implementation of these ideas using 
optically trapped molecules. We start by giving an es- 
timate of the infiuence of the tweezer potential which 
is used to guide the marker molecules. Since the marker 
and the crystal are made up of a similar type of molecules 
and spacings in our model are typically in the order of 
optical wavelengths the tweezer potential is also seen by 
molecules in the register. In a good approximation the 
potential of a tweezer focused at position (xq, ^o) is given 

by 



{x - xpf + {y- yof 
2aL 



(5.4) 



where Vq = muj^^cr^^ is the potential depth, <Jtw the 
width, and ujtw the resulting trapping frequency. To 
realize a controlled transition from a regime where the 
marker is only trapped by the tweezer to one where 
it is locked to a register molecule (see Fig. |3|, we re- 
quire that the width of the marker wave function in 
the tweezer potential is below the crystal lattice spac- 
ing, (fi/2mco'tw)"^^^ < ci' We expect little infiuence on the 
MDC if the dipole-dipole repulsion in the crystal exceeds 
the maximal force exerted on register molecules by the 
tweezer, 3y/eD/a'^ mu^^at^. These two conditions de- 
fine an allowed region for cjtw which gives the timescale 
for our adiabatic manipulations. 

For the example of LiCs molecules in an optical trap, 
induced dipole moments of about 3 Debye and a value 
of Td = 30 correspond to crystal lattice spacings of 
~ 630 nm. For perpendicular optical trapping frequen- 
cies ujj_ = 150 kHz such lattice spacings allow for min- 
imal values of b/a ~ 0.63 which can still be achieved 
with convenient laser wavelengths ~ 600 nm. In such a 
setup the strength of the dipole-dipole interaction in the 



crystal is 



(0) 
dd 



35 kHz such that reasonable ratios of 



Uo/Ures > 100 give absolute values for the local interac- 
tion of Uo ~ 10 kHz. The required deep optical trapping 
potentials can be realized with detunings of the order of 
100 GHz and Rabi frequencies of the order of 100 MHz. 
For red-detuned trapping fields this will lead to sponta- 
neous photon scattering rates of about Feff ~ 1 Hz which 
could be further reduced using blue detuned lasers. At 
the same time a tweezer that is focused to atw ^ l/J^m 
still allows for trapping frequencies uJtw between approxi- 
mately 0.6 kHz and > 6 kHz. This is compatible with the 
general time-scale of the effective spin-spin interaction 
and therefore adiabatic manipulations of the tweezer po- 
tential do not impose additional limitations on the speed 
of gate operations. 



We point out that the timescales in our model can 
be greatly enhanced for stronger perpendicular trapping 
frequencies. For example, by using electrostatic traps in 
on-chip setups trapping frequencies exceeding 1 MHz can 
in principle be achieved [27J. For such strong perpendic- 
ular confinement intermolecular spacings below 300 nm 
could be realized with reasonable ratios of b/a ~ 0.5, al- 
lowing for values of the local interaction strength easily 
up to 100 kHz while keeping Uo/Ures ^ 100. 

We finally comment on the readout of the quantum 
register. In contrast to atoms or ions the lack of closed 
optical cycling transitions in general prevents a non- 
destructive detection of molecular qubits. In previous 
proposals a state selective ionization of the molecules fol- 
lowed by detection of the charged molecule has been pro- 
posed as a readout scheme [22 . In our setting a similar 
method can be applied where in addition right before 
the readout stage a short time of free expansion of the 
crystal can be used to separate the molecules. Recently, 
a class of molecules with almost closed cycling transi- 
tions has been identified and used for laser cooling (47) . 
This would then also provide a non-destructive readout 
by first mapping the state of a register molecule onto 
such a 'closed-cycling' molecules using the gate schemes 
discussed above. 



VI. CONCLUSIONS & OUTLOOK 

In summary we have proposed and analyzed a new 
approach for implementing quantum gate operations be- 
tween molecular qubits embedded in a MDC. We have 
shown how arbitrary state dependent and time varying 
dipole moments for molecular spin qubits can be induced 
using microwave fields and how this can be exploited for 
the realization of enhanced phonon-mediated two qubit 
operations. Compared to previous schemes our proposal 
does not rely on a spectral or spatial resolution of in- 
dividual qubits and single site addressability is instead 
achieved with the help of marker molecules to modify 
the local phonon structure of the crystal. Therefore, this 
technique is in particular suited for molecules prepared 
in self- assembled MDCs, although it can be equally well 
applied to optical lattices or microtrap arrays. 

The general two qubit operation between marker and 
register molecules which we have analyzed here in more 
detail can serve as a basic building block for various scal- 
able quantum computing schemes with polar molecules. 
For an improved performance, different generalization 
can be envisioned, for example, the use of whole arrays of 
marker molecules to perform many operations in parallel 
or the multi-layer configurations to change from a 2D to 
a 3D architecture. Along these lines, also applications 
for digital quantum simulations [71 148] could be consid- 
ered, where the time evolution of arbitrary spin models 
is simulated by a stroboscopic application of two qubit 
gate operations. Here the fiexibility of moving marker 
molecules over several lattice sites between the opera- 
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tions can be used to adjust the effective range of the 
resulting spin-spin interactions. 
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